Self-consistent theory of phonon renormalization and electron-phonon coupling near a 

2D Kohn singularity 
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We show that the usual expression for evaluating electron-phonon coupling and the phonon 
linewidth in 2D metals with a cylindrical Fermi surface cannot be applied near the wave vector 
corresponding to the Kohn singularity. Instead, the Dyson equation for phonons has to be solved 
self-consistently. If a self-consistent procedure is properly followed, there is no divergency in either 
the coupling constant or the phonon linewidth near the offending wave vectors, in contrast to the 
standard expression. 
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First principles calculations of the phonon spectra and 
electron-phonon coupling in MgB 2 (see Ref. [1[ for a 
review) have determined that the interaction mainly re- 
sponsible for superconductivity in this material is cou- 
pling of small-g, high-energy optical phonons of a par- 
ticular symmetry with approximately parabolic nearly 
two-dimensional hole bands forming practically perfect 
circular cylinders occupying only a small fraction of the 
Brillouin zone. It was realized early enough Q that the 
case of an ideal 2D cylinder leads to a divergency in the 
calculated phonon linewidth at the 2D Kohn singularity, 
q = 2k p, which presents serious difficulties in calculating 
the electron-phonon coupling function. One option that 
was exploited was to use an analytical integration for the 
wave vectors comparable with, or smaller than 2kp 0, Q , 
and a numerical one for the larger vectors. 

It was also pointed out 0, |3| that this singularity gets 
stronger when the Fermi surface gets smaller, while the 
integrated electron-phonon coupling (for 2D parabolic 
bands) does not change. While a perfectly cylindrical 
Fermi surface is an idealized construction, deviations may 
be quite small, it seems, on the first glance, unphysical 
that all phonons with |q| = 2kp have infinite linewidth. 
Note that the problem is not specific for MgB2 : it oc- 
curs for any 2D material sporting Kohn singularities. In 
particular, the hypothetical hexagonal LiB, a subject of 
a substantial recent interest, has a bands that are even 
more 2D that those in MgB2 [|| and the described prob- 
lem is even more pronounced. 

This intuition is correct. In this paper we show that 
close to a Kohn anomaly standard formulas for calcu- 
lating electron-phonon interaction (EPI) become incor- 
rect, and new, self-consistent expressions replace them. 
These expressions have no singularities, and exhibit a 
much more natural, reasonably smooth, q-dependence of 
the phonon self-energy. 

To start with, we shall remind the readers the standard 
formalism. We first define the retarded phonon Green 
function 



D°"» (q,t) = -»0(t) 



uL(0) 



where [...] is a commutator and (...) denotes statisti- 
cal averaging. The displacement operator in the 
a— direction can be expressed via the phonon eigenvec- 
tors e q „ and frequencies squared u)^ u . 
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(i) 



For simplicity, a primitive lattice with a single kind 
of ions with a single mass M will be considered below. 
Also, atomic (Hartree) units will be used throughout the 
paper. In this case the "bare" phonon Green function 
has a form 



D (q,u) = 



1 



2Mw, 



1 



lu — luq + iS uj + ujq + id 



, (2) 



where luq is the bare phonon frequency, before account- 
ing for electron phonon coupling (screening by electrons) . 
Without losing generality, it can be assumed to be q in- 
dependent. Correspondingly, the full Green function is 



D(q,w) = 



1 



2MLU r 



(3) 

where w q is the renormalized (observable) frequency, and 
r q is damping (phonon linewidth) due to EPI 0. 

The Dyson equation reads 



^- 1 (q,^) = ^ - 1 (q,^)-n(q, W ), 



(4) 



where the polarization operator in the lowest approxima- 
tion (as usual, the Migdal theorem[7[ allows neglecting 
the vertex corrections) along the real frequency axis at 
T = has a form (a is the lattice constant in the plane) 



2 



n(q,w) = -2z y | ff °_ k+q | G (k + |, e + -)G (k - |, e - -)—— — . (5) 

Here k+q is the frare electron-ion scattering matrix element (the commonly used EPI matrix element differs in that 
the potential gradient is replaced by the derivatives with respect to the normal phonon coordinates) 



<?k,k +q = 4 / V k+q (r)W(r)^ k (r)d 2 r 

a J a 2 



(6) 



where 



G (k,e) = TTTT^T- — n 1 — T ( 7 ) 

e — eo(k) + iosign{k — kf) 

is the bare electron Green function, and £o(k) = (k\ + k 2 — k F )/2m. The renormalized phonon frequency and the 
phonon line width T q are determined by the pole of the phonon Green function -D(q, u>) or 

D _1 (q,o; q + ir q ) = L> ( 7 1 (q,tJ q + ir q )-n(q,cj q + ir q ) =0. 



This leads to 



and 



w q = w 2 + ^Ren(q, Wq +ir q )+r q , is; 



r q = -^imn ( q,. q + i r q ). (9) 

The next standard step, following Ref. Q, is to expand the polarization operator 

n (q ,„ + u, - -2S l< t ,J 2 eo(k + q 7/; ( y_ „ _ a do) 

to first order in frequency just above the real axis of the complex frequency u>. In this case 

= jj u ^fi) C 1 + °( w o/4)) (ii) 

and 

1 dlmll(q,w) 



papp _ 

q 2M dijj 



^|?El-9k,k +q | 2 ^o(k + q)]<5[£o(k)] (12) 

where the factor (nk — ™k+ q ) 5[eo(k + q) — £o(k) — lv] has been changed to —LuS[eo(k + q)]5[eo(k)]. According to Ref.Q 
"Except for extremely pathological energy bands, it is an excellent approximation". Unfortunately, MgB2 and some 
other recently discovered superconductors are examples where the energy bands are, in some aspects, pathological. 
Formally the expression (| 12(1 for 2D-system is divergent near a Kohn anomaly q — > 2kF, and we have r q > w q , i.e. 
phonons are not well defined quasiparticles. To describe electron-phonon interaction in these systems we have to 
calculate the polarization operator II(q, Z) for a complex frequency Z and solve Eqs. (|8|9p . 



COMPLEX POLARIZATION OPERATOR 



Let us consider a model with a cylindrical Fermi surface of radius kp whose electrons interact with an optical 
phonon with a bare frequency loq, and a momentum- independent matrix element g° (we also neglect possible warping 
of the Fermi-surface cylinder, cf. Ref. [loj ) . In this case the imaginary part of Eq. [TUl reads 

Imn(q, w + iS) = -2tt \g°\ 2 Y, {° M k + q)l - 9 Mk)]} 6 [e (k + q) - e (k) - w] 

k 
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ImII(q, u>) 



0| 2 ry/k 2 F +2mu 



2ixa 2 q 



Re 



dk 



max{fe J? ,ma)/g+«/2} J \ - (q/2k + mid/kq) 



ira 2 Q 



Rey/l - (Q - ft) 2 - Re^l - (Q + ft) 2 

(13) 



where Q = q/2kF, and 17 = lojqvF {vf is the Fermi velocity), 

To find the polarization operator for a complex frequency Z = lo + 17, we use the Hilbert transformation 



1 f 00 

n(q,Z) = -P/ __In 1 n(q ) J5). 

7T J_ 00 E-Z 



The result is 



H(q, Z) 



c dxdy -^/4e F )yQ» + x> : 



(14) 



where £p = k F /2m, x — k x /kp, y = k y /kp, and N(0) — m/2i:a 2 is the density of states at the Fermi level, per spin. 
The integration is performed over the range (x — Q) 2 + y 2 < 1 . The substitution x — Q = r cos ip, y — r sin <p leads to 



dxdy 



1 /-27T 

rdr I dip. 
Jo 



n( q , Z) = 



m \g" 



2tto? 



2-(l- 



4e F Q 2 



1-0- 



z 



4e F Q 



-{Z->-Z} 



(15) 



The branches of the square roots are chosen so as to get the correct behavior II oc q 2 Vp /lo 2 at large frequencies 
{lo = ReZ). For a 2D system for II(q, to) on the real frequency axis one can write 11, 13, 13, 14 1: 



ReU(q,uo + iS) = 



ImII(q, to + iS) 



m \g 



2ira 2 Q 

I 01 2 
m \g '\ 

' 2na 2 Q 



2Q-(Q- Q)Rcy/l-{Q-n)- 2 - (Q + Q)Rcy/l - (Q + fi)" 



Re-v/1 - (Q - ft) 2 - R ev /1 - (Q + ft) 2 



(16) 



First, we see that the imaginary part a finite for all values of the wavevector q and vanish inside the Landau-damping 
cone q < ui/vp (more exactly, at Q < lo/Aef — {lo/Aef) 2 )- It has two maxima: one is rather small, Imil(a;) ~ — 

m ^ 9 L^Zfff ; a t Q ~ l — id/4sp, while the other has an antiadiabatical behavior ImII(u;) ~ — "'^ 9 ^^J EF ^ 
occurs at a very low frequency Q ~ lo/Aef + (to/Asp) 2 . 

But if we expand the polarization operator at small frequencies we recover a standard result (see, e.g. Ref.0, 3) 



and 



n app (q,o; + iS) 



m \g 



12 r 



iid/Aep 



qVT~q 2 



(1 - Q) - y/1 - l/Q 2 6(Q - 1) 



(17) 



where the imaginary part diverges at q — > and q — > 2kp. The real part of Eq. 1161 practically coincides with the real 
part of Eq. [T7] except inside the Landau damping region. 
At id — > and finite q we get 



ImII a pp(q,w) 



™|<?°| 2 0(1 -Q) 



4e F 7ra 2 Q^J\ Z Qi 



(18) 



In the opposite limit ImII(q, to) = for q ^ mid/ ' \Jk 2 F — mid. For the real part one can set ui = in Eq. 1171 



Ren(q, 0) = -2 \g°\ 2 ^ 0(|k| - k F )9(k F - |k + q|) 



e (k + q) -eo(k)' 
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In this case 



Ren(q, 0) ~ - 
In the opposite limit 



m \g 



1 - 0(q - 2k F )^l-{2k F /qf 



m \g 



1 - 



-l)v/l-Q- 2 



(19) 



ReII(q-» 0,lu) 



m \g 



2na 2 



qk F 
raw 



The momentum dependence of the absolute values of the imaginary (the upper panel) and the real parts (the 
lower panel) of Eqs. [16] (solid lines) and [17] (short-dash lines) at oj = ojq = 90 meV as the functions of the reduced 
wavevector q/ksz (ksz = t/o is the Brillouin zone vector, or the radius of the Wigner-Seitz cylinder) are shown in 
Fig. 1. Fermi vectors k F /kBZ = 0.075, k F /kBZ = 0.1,/cf/&bz = 0.15, and kpjksz — 0.2 correspond to e F =0.15 eV, 
0.27 eV, 0.60 eV, and 1.07 eV, respectively. 

Along the imaginary (Matsubara) axis the polarization operator has the following form 



n M (q,iWn) 



m \ g 



V2Q 2 



(20) 



where uo n — 2imT. T is temperature, n=0, ±1, ±2, ... ± oo. IlM(q, i^n — ► 0) coincides with the Eq. 1191 

! — 



PHONON RENORMALIZATION IN 2D 
SYSTEMS 



First let us consider the approximate polarization op- 
erator from Eq. [T7] Then 



0J, 



app 



w 0A /l-2C 



1 - 6{q - 2k F )^l - (2k F /qf 



and according to Eqs. [T21 and [T7l 



(21) 



j^app 

q 



m \g 



6(2k F ~q) (to 2 9(1 -Q) 



2-KO?MV F q y/i - (q/2k F ) 2 4e F Q^l-Q 2 ' 

n (22) 

where we have introduced, following Ref. [1 51 ] . an aux- 
iliary coupling constant C = -^(0) 3° /Mco 2 (some 
authors use another dimcnsionlcss constant Ao = 2£). 
il q = o>q PP + iT^' p gives the renormalized frequency and 
the damping (see Eqs. Ill|12|) . For q < 2k F 

1 



T^app 

q 



(23) 



This expression diverges in the limits Q — > and Q — > 1. 

Turning now to the exact Eq. I15[ we observe that in 
the small q limit the polarization operator becomes real: 



n(q^o, w ) 

Solving the equation 



m \g" 



/ 4e F 
27ra 2 \ ui 



Q 2 



(24) 



£> _1 (q ^0, w) = w 2 - + C^o 



= 



one gets 




16Cwg4Q 



We choose the branch of the complex square root that 
gives Wq^o - > for £ — > 0. This means that up to 
0(Q 2 ) there is no damping for the phonon. The phonon 
spectral function F(u>) = — -InxD(q — » 0, w) shows a 
narrow peak at this frequency and Fq^o/^q^o vanishes. 
The situation at Q — > 1 is similar. In the lowest order in 
lu/e f we have 



Tl(q^>2kp,u) = 



m \g 



e F 



This leads to 



and 



^2k F 
^2k F 

^2k F 



04 

^2k F 



•2C, 

0J2k F 

8e F 



( woVT^2C 



^2k f 



1-2C 



8e? 



(25) 



This ratio remains finite in the limit Q — > 1, although the 
approximate expression of Eq. 1231 diverges for any system 
with a cylindrical Fermi surface . The main point is that 
in both cases the well known popular formula 



k 



(26) 



k 2 l 



Af(27ra) 



dk 



|v(k) x v(k + q) 
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FIG. 1: (color online) The imaginary and real parts of the normalized polarization operator 2% |II(q, uj = u>o) \ a 2 /m \g°\ as a 
function of the reduced wavevector q/ksz, for different fillings kp/ksz- Solid lines represent the exact results, and the dashed 
lines the approximate solution (Eq. I17jl . The four different sets correspond, from left to right, to four increasing values of 

kF /kBZ- 



is not valid near the Kohn singularity. 

The results for the linewidth T q and the renormal- 
ized phonon frequency w q obtained by using the approx- 
imate polarization operator as functions of the reduced 
wavevector Q = q/2kp are shown in Fig. 2 by red lines. 
Parameters are following : the bare phonon frequency 
ljq = 90 meV , the bare constant of EPI ( = 1/4. The 
ratio kp/kBz is equal to 0.17. It corresponds to Ef — 0.2 
eV. 

The approximate result agrees with that of Ref.Q 
(their Fig. 4). The exact result (black lines) has been 
obtained by a numerical solution of Eqs. I8|9I using the 
polarization operator from the Eq. 1151 The latter, in 
contrast to the approximate expression, shows two shoul- 



ders in the wavevector dependence of the renormalized 
frequency w q (see also |ReII(q, ujq)\ in the bottom panel 
of Fig. 2). The first one corresponds to the maximum of 
|ImII(q, w)| (or T q ) and the second one to vanishing of 
these values. 



ELECTRON SELF-ENERGY 

The electron self-energy is expressed via the electron 
and phonon Green functions: 



E(p) = ^.g / G(p - k)D(k)T(p,p - k, k) 



d D+1 k 



6 



> 



5=0.25, 




^ co o =90 meV / 


exact solution 


approximate solution 








where p = {p,s}. It was shown by Migdal 0| that 
one can neglect the vertex corrections T(p,p — k,k) ~ 
g°(l + y/m/M) and that the function G(e, k) differs from 
the bare electron Green function (Eq. [7]) only in the nar- 
row interval of momenta \k — kp\ < uj p h/VF and frequen- 
cies |cj| < ij^ph- Thus, the full electron Green function 
G(p) can be substituted by the corresponding function 
for noninteracting electrons (Eq. [7]). Using the Eq. [3l 
the electron self-energy for T = becomes (see, e.g., 
Refs.0, 



q/2k F 



FIG. 2: (color online) The linewidth T q and the renormalized 
phonon frequency u q obtained by using the approximation of 
Eq. [T7] and the exact expression Eq. fK] for the following 
parameters : the bare phonon frequency u>o = 90 meV , the 
bare constant of EPI £ = 1/4. The filling corresponds to 
k F /k B z = 0.17 and e F = 0.2 eV. 



S (k,.)=^ (£kW Jfi/ 



0(0 



This is a trivial generalization of the standard expressions onto the finite phonon linewidth case. Let us average the 
self energy over the Fermi surface 

E H = ^yE^)s(k,^) 



2^ 2^ r5(e k )(5( £k+q )^ ^ In ^ ; - + 2i 



N(0)M 



k q 



r q + K + ^ 



tan 



tan 



The limit A = — lim^^o ReE(o;)/oj is nothing but the standard electron-phonon coupling constant 



MiV(O) 



k q 



m r q + < 



(27) 



where we introduced the phase space function (sometimes called "nesting function" ) 



A. 



,(°) = ]v|5)E^)^k +q )- 



(28) 



For a 2D cylindrical Fermi surface we have 



^q(O) 



(1-0) 



(29) 
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which diverges at Q — > and Q — > 1. Following Ref. [8| we can introduce the "mode A" via the expression A = Y^ q ^q- 
Then 



A q = iV q (0)^ 



1 



m r q + c q - 



In the weak-damping approximation for we recover the standard formula 

I °| 2 

A£*=JV q (0)!£J-/a& 

but for T q 3> w q the contribution of strongly damped phonons to total A q is suppressed. 

The result (f37|) we can get also if we introduce, according to Eq. [31 a generalized Eliashberg function 



al(uj)F(uj) 



*(e k )|s°| *(ek+q)l 



27rAfiV(0) 



-E 



2ttM 



E 

k,q 

^ q (o)|g°| 2 i 

2w q 2 



2ui c 



(W - CJq) 2 + T q ( W + CJq) 2 + r 



(30) 



(31) 



The second term in this expression cancels out the nonphysical behavior at low and high frequencies. Otherwise Ar 
would have be divergent. Eqs. 1271311 are general and valid for not only for the 2D systems, where phase space factor 
(Eq. is divergent. 



For T q <C w q we have 



— y 



5(e k )|ff°r<J(e k+q ) 



N(0)M 



k,q 



2u r , 



1 y°-pp 



where in the last equality we have used the approximate Eq. [T2l Note that the Eq. [32] is a consequence of that 
fact that the damping r q pp , according to Eq. [T^l can be expressed via the "nesting function" (Eq. |2"5]) and both are 
determined by the same function ImII'(g,0) = dhxJI(q,u))/duj\u=o- In a general case these functions can be different. 
This result without using the pole approximation for the phonon Green function can be trivially obtained in the 
Matsubara formalism. In this case it one does not need to solve the Dyson equation. In the lowest order in coupling 
for T = the self-energy has a form 



ilu — £k' — S(icJ, k') ' 



where 



%„(!w,k,k') = (1/tt) / dnimD(Q + i8, k, k') \(iw - ft) 1 - (iw + ft) 1 
We can also average the self-energy over the Fermi surface S(ie) = Y^ k i5(ek)S(ie, k)/iV(0). 



(33) 



1 f'OC J p 

^) = ^ EE ^I^I^W J^—D^iu) J 



k q 



-oo i(e-w) -£ 



The integral J_ l30 — — ?27rsign(e — u>) allows to calculate the physical coupling constant A 

k q 



A 



die 



]™^EE^)l.9°l 2 ^ k+q ) / ^D^i^ 



S(e - to) 



«jo)EE <*(<«<> l»T '<«*JD -«■• *» - -4) EE j m l*"!"'^^^ 



k q 



k q 



D M L o (i0)~U M (q,i0) 



r 



On the Matsubara axes, for 2D system, according to Eq. Hm(<1 < ^kp, iO) = —2(uiq is a constant (see Eqs 12171 and 
[20l the phase space factor vanishes for q > 2k p and 
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US). Using Eqs. I2"l33l we get D^ (iO) = -Mwg and 



A = C/(1-2C). 



(34) 



In a 3D case Iljvf (q, iO) is a rather complicated function 
of q and the Eq. [M] is only an approximation (as was 
probably firstly mentioned by Frohlich 15|). The physi- 



cal meanings of the coupling constants A and £ that they 
are measures of the renormalization of the phonon fre- 
quency from wo to w q (c/. a discussion for 3D systems in 

Ref.Q )• 

Turning back to the 2D case, according to Eq. [57] we 
can neglect all divergent contributions near Q = and 
Q = 1. Elsewhere we can use Eq. [30l 

One should keep in mind that the conventional cou- 
pling constant isA = £/(l — 2£). This parameter deter- 
mines electronic properties (Fermi velocities, T c , etc.). 
The other parameter, £ = A/(l + 2A) < 1/2, defines the 
observable phonon frequency, w q = u>o\/l — 2£. 



CONCLUSIONS 

First of all, the standard well-known expression 
a 2 (u))F(uj) 
1 



AT(0)M 



^7 *(«-«q), 



k q 



is valid only in the lowest order in the "bare" phonon 
linewidth, r q , which is not an acceptable approximation 
in case of strong Kohn singularities, and particularly for 
a cylindrical Fermi surface. T q in this approximation is 
not the actual phonon line width; as opposed to r q , de- 
termined by the oversimplified Eq. I12( the real phonon 
linewidth does not diverge even for an ideally cylindrical 
Fermi surface. Second, the renormalized A for a cylindri- 
cal Fermi surface and parabolic bands does not depend 
on filling. 
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